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It is argued by the author that the canonical form of the quark energy-momentum tensor with a 
partial derivative instead of the covariant derivative is the correct definition for the quark momen- 
tum and angular momentum fraction of the nucleon in covariant quantization. Although it is not 
manifestly gauge invariant, its matrix elements in the nucleon will be non-vanishing and are gauge 
invariant. We test this idea in the path-integral quantization by calculating correlation functions 
on the lattice with a gauge-invariant nucleon interpolation field and replacing the gauge link in the 
quark lattice momentum operator with unity, which corresponds to the partial derivative in the 
continuum. We find that the ratios of three-point to two-point functions are zero within errors for 
both the u and d quarks, contrary to the case without setting the gauge links to unity. 



There is an ongoing controversy in QCD over the is- 
sue of how to separate total angular momentum between 
the quark and glue sectors and, moreover, whether there 
exists a gauge-invariant decomposition of the glue angu- 
lar momentum into spin and orbital angular momentum 
components as in the quark case. While several versions 
of gauge-invariant separation of the glue spin and orbital 
angular momentum have been proposed [TJSj, X. Ji ar- 
gued that such a separation cannot be be achieved [5] . On 
the other hand, E. Leader asserts [6] that, if the momen- 
tum operators are to be the generators of spatial trans- 
lations, then the gauge-invariant Bcllinfante momentum 
operators are not correct and that rather it is the canon- 
ical form that is correct. 

Since this approach is based on covariant quantiza- 
tion, we shall use the path-integral quantization to check 
this conclusion. In the Monte Carlo simulation of QCD 
on the lattice, the number of degrees of freedom is fi- 
nite, and one does not need to fix the gauge. In this 
case, it is usually presumed that gauge- variant observ- 
ables vanish according to Elitzur's theorem. However, 
E. Leader suggests [7] that, despite the fact that the op- 
erator is not manifestly gauge invariant, there could be 
a gauge- invariant part of the matrix element. To check 
this possibility, we perform a lattice gauge calculation of 
the quark energy-momentum tensor operator. To calcu- 
late the quark momentum fraction in the nucleon on the 
lattice, one typically considers the ratio of the 3-point to 
2-point function 

1 r Ti:YGNON{tf,t,to\p l ) tf _ t> i,t_t »i / y ,y 
Pi TrTG NN {tf,t ;pi) 

where T is the projector onto the unpolarized nucleon 
state. The two-point and three-point functions for the 
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nucleon are 

^e-^-^ix^m^xhM), (2) 

X 

J- e-^^\x a {x f ,t f )0{x,t)^(xM), (3) 

Xf,X 

where \ is the nucleon interpolating field and to A/ is 
the source/sink time of the nucleon interpolation field. 
The lattice operator O for the momentum that has been 
calculated in the literature jS] is 

Tu{x) = ^Lft^UiixW + ...} (4) 

with Ui(x) being the gauge link. This corresponds to 
the continuum tensor operator for the quark momcn- 

turn Pi = ^ipijiDi +jiDi)^p which is of the Bcllinfante 
form. E. Leader argues [6j that the correct momentum 
operator should be the canonical form, with the covari- 
ant derivative in the Bellinfante form replaced with the 
partial derivative. This corresponds to setting the gauge 
link U to unity in the lattice operator Ta(x) in Eq. Q. 

We performed a lattice calculation on a 2 + 1 + 1-flavor 
24 3 x64 HISQ lattice [S] (HYP-smeared) with the Wilson- 
clover valence fermions. The pion mass is 310 MeV with 
lattice spacing around 0.12 fm, and 127 configurations 
(each with 4 sources) were used in the calculation. The 
nucleon source is shifted to to = and sink timeslice is at 
tf = 12. The ratio in Eq. |l]) for the connected insertion 
is plotted in Fig. [I] as a function of the operator insertion 
time t for the u quark in the upper panel and the d quark 
in the lower panel. The results with the Bcllinfante form 
are shown as blue circles; they are obviously non-zero. 
The quark momentum fraction (x) is usually extracted 
from a plateau in a window away from the source and the 
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sink. On the other hand, when the gauge link U is set 
to unity, which corresponds to a partial derivative, the 
results (red squares) are zero within errors for practically 
all insertion times t. We conclude from this calculation 
that the canonical form of the quark energy-momentum 
tensor with partial derivatives does not lead to non-zero 
matrix elements in the nucleon. 
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FIG. 1. The ratio of the three-point to two-point functions in 
Eq. |l]) is plotted as a function of the operator insertion time 
t for the u/d quark in the upper/lower panel. The blue circles 
are results from the momentum operator with the covariant 
derivative and the red squares are those corresponding to the 
partial derivative. 
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